Scattering length of the helium atom - helium dimer collision 
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We present our recent results on the scattering length of '*He-'*He2 collisions. These 
investigations are based on the hard-core version of the Faddeev differential equations. As 
compared to our previous calculations of the same quantity, a much more refined grid is 
employed, providing an improvement of about 10%. Our results are compared with other 
ab initio, and with model calculations. 
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I. INTRODUCTION 

Weakly bound small ^He clusters attracted considerable attention in recent years, in particular 
because of the booming interest in Bose-Einstein condensation of ultra-cold gases 

Experimentally, helium dimers have been observed in 1993 by Luo at al. [jsl], and in 1994 by 
SchoUkopf and Toennies |4]. In the latter investigation the existence of helium trimers has also 
been demonstrated. Later on, Grisenti et al. measured a bond length of 52 ± 4 A for ^He2, 
which indicates that this dimer is the largest known diatomic molecular ground state. Based on 
this measurement they estimated a scattering length of 104^jg A and a dimer energy of l.l+jj;3mK 
1 5]. Further investigations concerning helium trimers and tetramers have been reported in Refs. 
1 6,0], but with no results on size and binding energies. 

Many theoretical calculations of these systems were performed for various interatomic poten- 
tials [8, 9, 10, 11]. Variational, hyperspherical and Faddeev-type techniques have been employed 
in this context (see, e.g., lll2ll - ll28l] and references therein). For the potentials given in ||9|,ll0|] it 



turned out that the Helium trimer has two bound states of total angular momentum zero: a ground 
state of about 126 mK and an excited state of about 2.28 mK. The latter was shown to be of Efimov 
nature I.13, .15.. .211. In particular, it was demonstrated in |21] how the Efimov states emerge from 
the virtual ones when decreasing the strength of the interaction. High accuracy has been achieved 
in all these calculations. 

While the number of papers devoted to the "^Hes bound-state problem is rather large, the number 
of scattering results is still very limited. Phase shifts of ^He-^He2 elastic scattering at ultra- 
low energies have been calculated for the first time in iIitL [Till below and above the three-body 
threshold. An extension and numerical improvement of these calculations was published in 1E4I1 . 
To the best of our knowledge, the only alternative ab initio calculation of phase shifts below the 
three-body threshold was performed in ['28']. As shown in |2^ sO], a zero-range model formulated 



in field theoretical terms is able to simulate the scattering situation. 

Though being an ideal quantum mechanical problem, involving three neutral bosons without 
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complications due to spin, isospin or Coulomb forces, the exact treatment of the He triatomic 
system is numerically quite demanding at the scattering threshold. Due to the low energy of the 
Helium dimer, a very large domain in configuration space, with a characteristic size of hundreds 
of Angstroems, has to be considered. As a consequence, the accuracy achieved in 1 19, 24] for the 
scattering length appeared somewhat limited. To overcome this limitation, we have enlarged in the 
present investigation the cut-off radius Pmax from 600 to 900 A and employed much more refined 
grids. 



II. FORMALISM 



Besides the complications related to the large domain in configuration space, the other source 
of complications is the strong repulsion of the He-He interaction at short distances. This prob- 
lem, however, was and is overcome in our previous and present investigations by employing the 
rigorous hard-core version of the Faddeev differential equations developed in I31 ..,32.1. 

Let us recall the main aspects of the corresponding formalism (for details see [191 12411). In what 
follows we restrict ourselves to a total angular momentum L = 0. In this case one has to solve the 
two-dimensional integro-differential Faddeev equations 
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0, X < c . 



(1) 



Here, x,y stand for the standard Jacobi variables and c for the core range. The angular momentum 
I corresponds to a dimer subsystem and a complementary atom; for an S-wave three-boson state, 
I is even (I = 0,2,4, . . .). V{x) is the He-He central potential acting outside the core domain. The 
partial wave function is related to the Faddeev components 4>/(x, by 



'i'i{x,y)=^i{x,y) + Y, / dr]hnix,y,r])%ix,y') 
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and 1 < 7] < 1. The explicit form of the function hui can be found in Refs. 11331 134|] . 
The functions 4>/(x,j) satisfy the boundary conditions 
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Moreover, in the hard-core model they are required to satisfy the condition 

r+l 

j ^ drihii^{c,y,r])^,^{x,y')=0. 



(3) 



(4) 



This guarantees the wave function to be zero not only at the core boundary x = c but also 

inside the core domains. 

The asymptotic boundary condition for the partial wave Faddeev components of the two- 
fragment scattering states reads, as p ^ oo and/or y ^ <>°, 



^i{x,y\p) 



dio\irdix)^sm{py)+exp{ipy) ao{p)+o{y 
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(5) 



3 



Here, y^di^) is the dimer wave function, E stands for the scattering energy given hy E = + p 
with £j the dimer energy, and p for the relative momentum conjugate to the variable y. The 

variables p = Wx^+y^ and G = arctan- are the hyperradius and hyperangle, respectively. The 

X 

coefficient ao{p) is nothing but the elastic scattering amplitude, while the functions A/(0) provide 
us, at £■ > 0, with the corresponding partial-wave Faddeev breakup amplitudes. The ^He - '*He2 
scattering length isc is given by 

= i!2M. (6) 

2 p 

Surely we only deal with a finite number of equations dU)-©, assuming / < l^ax, where /max 
is a certain fixed even number. As in [19, 24] we use a finite-difference approximation of the 
boundary-value problem ([T])-© in the polar coordinates p and 6. The grids are chosen such that 
the points of intersection of the arcs p = pi, i = 1,2, ... ,Np and the rays G = Oj, j = 1,2, ... ,Nq 
with the core boundary x = c constitute the knots. The value of the core radius is chosen to be 
c = 1 A by the argument given in [24]. We also follow the same method for choosing the grid radii 
pi (and, thus, the grid hyperangles Gj) as described in 1.1 9. .24.1 . 



in. RESULTS 



Our calculations are based on the semi-empirical HFD-B and LM2M2 potentials by 
Aziz and co-workers, and the more recent, purely theoretically derived TTY [ 10] potential by 
Tang, Toennies and Yiu. For the explicit form of these polarization potentials we refer to the 
Appendix of Ref. 112411 . As in our previous calculations we choose h^/m = 12.12 KA^, where 
m stands for the mass of the ^He atom. The ^He dimer binding energies and ^He-^He scattering 
lengths obtained with the HFD-B, LM2M2, and TTY potentials are shown in TableHl Note that the 
inverse of the wave number x^^) = a/IcTI lies rather close to the corresponding scattering length. 



TABLE I: Dimer energy Ed, wave length l/x^^', and ^He— '^He scattering length llJ for the potentials used, 
as compared to the experimatal values of Ref. [5]. 



£j (mK) iil^ (A) 


Potential 


Sd (mK) 


l/x(2) (A) 


(A) 




LM2M2 


-1.30348 


96.43 


100.23 


Exp. [5] l.llj^i 104+^jj 


TTY 


-1.30962 


96.20 


100.01 




HFD-B 


-1.68541 


84.80 


88.50 



Unlike the trimer binding energies, the ^He-^He2 scattering length is much more sensitive 
to the grid parameters. To investigate this sensitivity, we take increasing values of the cut-off 
hyperradius pmax, and simultaneously increase the dimension of the grid N = Nq = Np. Surely, in 
such an analysis we can restrict ourselves to /max = 0. The results obtained for the TTY potential 
are given in Tablelllland Fig. [T] Inspection of this figure shows that, when increasing the dimension 
N of the grid, convergence of the '^He-'^He2 scattering length isc is essentially achieved, however, 
with different limiting values of E^c for different choices of Pmax- This concerns, in particular, the 
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TABLE II: The '^He-'^Hei scattering length ^sc (A) for 

^max — in case of the TTY potential as a function 

of the grid parameters pmax and N = Np =Ne- 



Pmax 


1005 


1505 


2005 


2505 


3005 


3505 


600 


162.33 


159.80 


158.91 


158.61 


158.31 




700 


164.13 


159.99 


158.57 


157.99 


157.65 


157.48 


800 


167.15 


160.98 


158.90 


158.03 


157.46 




900 


171.19 


162.52 


159.66 


158.40 


157.66 





TABLE III: The '^He-'^He2 scattering length £,c (A) obtained for a grid with Np = iVe=2005 and p^^^=700 A. 



Potential 


^max 


This work 


[24] 


[23] 


[28J 


[35J [30] 







158.2 


168 








LM2M2 


2 


122.9 


134 










4 


118.7 


131 


126 


115.4 


114.25 113.1 







158.6 


168 








TTY 


2 


123.2 


134 










4 


118.9 


131 




115.8 


114.5 







159.6 


168 








HFD-B 


2 


128.4 


138 










4 


124.7 


135 




121.9 


120.2 



transition from Pmax = 600 A to Pmax = VOO A, while the transition to 800 A or even 900 A has 
practically no effect. 

Bearing this in mind, we feel justified to choose Pmax = 

700 A when going over from Uax = 
to /max = 2 and 4. The corresponding results are presented in Table |III| There we also show 
our previous results [12411 where, due to lack of computer facilities, we had to restrict ourselves to 
Pmax = 460 A and A'^ = 605. We see that an improvement of about 10% is achieved in the present 
calculations, as indicated already by the trends in Fig. [H 

Table HiH also contains the fairly recent results by Blume and Greene iEbI i and Roudnev lEsll . 
The treatment of [23] is based on a combination of the IVIonte Carlo method and the hyperspherical 
adiabatic approach. The one of Ref. lEsll employs the three-dimensional Faddeev differential 
equations in the total angular momentum representation. Our results agree rather well with these 
alternative calculations. 

This gives already a good hint on the quality of our present investigations. A direct confirmation 
is obtained by extrapolating the curves in Fig. [T] According to this figure, convergence of isc as 
a function of is essentially, but not fully, achieved. A certain improvement, thus, is still to be 
expected when going to higher A''. In order to estimate this effect we approximate the curves of 
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FIG. 1 : The '^He-'^He2 scattering length ^sc for l„ 
grid dimension N = Np =N0. 



in case of the TTY potential as a function of the 



Fig. [Uby a function of the form 

UN) = a + j^^. (7) 

Clearly, £sc{°°) = GC. The constants a, /3, and 7 are fixed by the values of £sc at A^^ = 1005, 2005, 
and 3005. In this way we get the corresponding optimal scattering lengths £sc(°°) = 157.5, 156.4, 
155.4, and 154.8 A for pmax = 600, 700, 800, and 900 A, respectively. Comparing with Table 
shows that the differences between these asymptotic values and the ones for = 3005 lie between 
1 to 3 A. 

For /max = 4, pmax = 700 A and the LM2M2 potential the scattering length has been calculated 
for A'^ = 1005, 1505, and 2005. Employing again the extrapolation formula Q with a, /3, 7 being 
chosen according to these values, we find £sc(°°) = 1 17.0 A. The difference between the scattering 
length obtained for A'^ = 2005 and the extrapolated value, hence, is 1 .7 A. A direct calculation for 
higher A'^ should lead to a modification rather close to this result. Following this argumentation, 
we conclude that the true value of for the LM2M2 and TTY potentials lies between 115 and 
116 A. 

For completeness we mention that besides the above ab initio calculations there are also model 
calculations, the results of which are given in the last two columns of Table |llll The calcula- 
tions of llisll are based on employing a Yamaguchi potential that leads to an easily solvable one- 
dimensional integral equation in momentum space. The approach of llioll (see also i29ll and ref- 
erences therein) represents intrinsically a zero-range model with a cut-off introduced to make the 
resulting one-dimensional Skomyakov-Ter-Martirosian equation [36] well defined. The cut-off 
parameter in |29, 30] as well as the range parameter of the Yamaguchi potential in [35] are ad- 
justed to the three-body binding energy obtained in the ab initio calculations. In other words, these 
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approaches are characterized by a remarkable simplicity, but rely essentially on results of the ab 
initio three-body calculations. 
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